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Two known 2-dim SUSY quantum mechanical constructions - the direct gen- 
erahzation of SUSY with first-order supercharges and Higher order SUSY 
with second order supercharges - are combined for a class of 2-dim quantum 
models, which are not amenable to separation of variables. The appropriate 
classical limit of quantum systems allows us to construct SUSY-extensions 
of original classical scalar Hamiltonians. Special emphasis is placed on the 
symmetry properties of the models thus obtained - the explicit expressions of 
quantum symmetry operators and of classical integrals of motion are given 
for all (scalar and matrix) components of SUSY-extensions. Using Grass- 
manian variables, the symmetry operators and classical integrals of motion 
are written in a unique form for the whole Superhamiltonian. The links of 
the approach to the classical Hamilton- Jacobi method for related "flipped" 
potentials are established. 
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1. Introduction 



Supersymmetric Quantum Mechanics (SUSY QM) Ij is a new framework for analyzing 

non-relativistic quantum problems. In particular, it helps to investigate the spectral prop- 
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erties of different quantum models as well as to generate new systems with given spectral 
characteristics (quantum design). 

Much less attention has been paid to SUSY QM as a tool to study individual hidden 
symmetries of the superpartner Hamiltonians. This problem is reasonable for considering 
either one-dimensional quantum systems with internal degrees of freedom (with matrix po- 
tentials) or systems of higher spatial dimensionality. Only for these classes of systems may 
the symmetry operators, which are in involution with the Hamiltonian (and are independent 
of it), exist. Thus, the quantum models in one- dimensional SUSY QM with matrix poten- 
tials and the higher-dimensional (in particular, two-dimensional) models with scalar and/or 
matrix potentials are extremely attractive. 

Regarding SUSY QM systems with matrix potentials, we refer to papers The SUSY 
QM systems with an arbitrary {d > 1) dimensionality of space'^ were constructed and in- 
vestigated in [S]. Such models include both the scalar and matrix components of the 
Superhamiltonian. The latter are interesting either for a description of interacting non- 
relativistic particles with spin [Hj, [Zj, or for developing supersymmetric quantum field the- 
ory on a spatial lattice 0. The appearance both of scalar and matrix Hamiltonians in a 
unique Superhamiltonian provides an opportunity to consider (starting from a given scalar 
Schrodinger operator) SUSY extensions that correspond to the systems with internal degrees 
of freedom. Some SUSY extensions of this type were considered in [HI (Calogero-like models 
of particles on a line) and in (Coulomb potential in d dimensions). 

Among multidimensional SUSY QM models, those most developed are the two- 
dimensional ones. Namely, precisely for these systems second-order supercharges were used to 
build the higher-order deformation of SUSY algebra JT], ^21; CSl- In the framework of this 
HSUSY QM generalization of conventional Witten's superalgebra'^ one can avoid the appear- 
ance of any matrix components of the Superhamiltonian, so that two scalar two-dimensional 
Schrodinger operators are intertwined by second-order supercharges. As a by-product of this 
construction, each of the intertwined Hamiltonians obeys the hidden symmetry: the differen- 
tial operators of fourth order in derivatives exist, which are not reducible to the Hamiltonian 
and commute with the Hamiltonian, ^1], ^2], ^Ej, jlB]- In the two-dimensional context this 
■^Instead of dimensionality of the space d can be also interpreted as a number of particles, for example in 

Calogero-like models of interacting N particles on a line [3] . 

''One-dimensional HSUSY QM was investigated in detail in \TE[. 
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means the complete integrability of the system. 

Another direction in which to investigate SUSY QM models involves their connections 
with classical counterparts. Initially, the arbitrary-space-dimensionality generalization of 
SUSY QM, mentioned above, was obtained by canonical quantization of a suitably chosen 
multi-dimensional system of Classical Mechanics. Then the quasi-classical limit of some su- 
persymmetrical quantum models investigated afforded new insight into the properties of the 
classical models obtained. In the one-dimensional case, this limit led to a new SWKB quan- 
tization rule ^Tj, which turned out to be more useful than the standard WKB rule. In the 
two-dimensional case, the quasi-classical limit provided an alternative effective method JH] 
for the construction of integrable systems in Classical Mechanics, which essentially enlarges 
the list of such models. The SUSY QM approach also provides new interesting links with 
the well-known Hamilton- Jacobi equation in Classical Mechanics 

In the present paper we shall combine both known two-dimensional SUSY QM construc- 
tions - the direct generalization of SUSY with first-order supercharges and HSUSY with 
second-order supercharges - in order to investigate the symmetry properties of the models 
obtained, both at quantum and classical level. The paper is organized as follows. In Sec- 
tion 2, known results about two-dimensional SUSY QM models and their connections with 
classical models, necessary for the original part of the paper, are briefly summarized. In 
Section 3, the particular case of two-dimensional models for second-order supercharges with 
intermediate twists are studied: the particular models of generalized Morse and Poschl- Teller 
potentials are presented within this common framework. In Section 4, the integrability of 
these models is extended onto their matrix - both quantum and classical - superpartners. In 
Section 5, the links with classical Hamilton- Jacobi equation are considered. 

2. Two-dimensional SUSY Quantum Mechanics 

2.1. The representation of SUSY algebra with first-order super- 
charges 

In the two-dimensional case x = (xi,X2), the direct (i.e. of first order in derivatives) 
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generalization of SUSY QM satisfies j3] the conventional Witten's [20] SUSY algebra 



{g+, Q-} = hH- {g+, g+} = {q-,q-} = o; [g±, ^] = o. 



by the 4x4 matrix operators: 
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where two scalar Hamiltonians H^^\H^'^^ and one 2x2 matrix Hamiltonian of the 



Schrodinger type take on a quasi-factorized form: 

H(') = q+q- = -h'df + = -h'df + (daix))' - hdfxix); df ^ df + 



(3) 
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fi'S^kdf + 5^J{dixix)Y - hdfxix)) + 2hd,dkx{x 



the components of the supercharges being of first order in derivatives: 



qf = ^hdi + (dixix)); sf 



(4) 



where di = d/dxi and summation over repeated indices is assumed. The Planck constant 
was restored in formulas (01-0, and the {normalizable or unnormalizahle) zero-energy 
wave functions of the scalar Hamiltonians are now written as exp (±x/^). 

The quasi-factorization in (jH)) ensures that the last equation in holds, and leads to 
the following intertwining relations for the components of the Superhamiltonian (P)): 



9. 



r(l) 



(1), 



r(i) 

ki ■ 



(5) 



These relations play the main role in the SUSY QM approach and, in particular, they lead 
to the connections between the spectrum of the matrix Hamiltonian and the spectra of two 
scalar ones. We remark that i^'-^^ and if'-^^ are not intertwined with each other and are not 
(in general) isospectral since q'^ ■ =0. 
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2.2. Second-order supercharges 



An alternative opportunity to include two-dimensional scalar Hamiltonians in the SUSY 
QM framework is based on the supercharges of second order in derivatives ^T], jT2j, P^ : 

= {Q')^ = gik{x)tedidk + Ci{x)hdi + B{x) (6) 

where gik^ Ci, B are arbitrary real functions. In this case, two scalar Hamiltonians H^'^\ H^^^ 
are intertwined directly without any (matrix) intermediary: 

H(^\x)Q+ = Q+H^^^x); H^''\x)Q- = Q-H^''\x). (7) 

Although no method to find the general solution of the intertwining relations ((Tj) has been 
proposed, a certain number of models obeying these relations have been found for the cases 
of hyperbolic (Lorentz) and degenerate metric gik- One very important specific property of all 
these models, which follows from the intertwining relations (jTj), is their integrability. Indeed, 
both Hamiltonians possess the symmetry operators R^'^\ R^^^ of fourth order in derivatives 

ijW] = 0; [R^^\ #(°)] = 0; i?^ = Q'Q^; R^^^ = Q+Q-, (8) 

which are not, in general^, polynomials of H^'^\ H^^\ 

For the case of Lorentz (hyperbolic) metric g^k = diag{l, —1) [llj - the intertwining 
relations (jZj) can be rewritten in a reduced form: 

d4C^F) = -d+{C+F), (9) 

where x± = {xi ± X2)/V^; functions Ci^2 were found to satisfy C± = Ci ^ C2 = C±{\/2x±), 
and F{x) = Fi{2xi) + ^2(2x2). Thus, the potentials and the supercharges are 

expressed in terms of the functions C±{^/2x±) and Fi(2xi), ^2(2x2): 

V'^^\V^''> = Th(c'4V2x+) + C'_{V2x.)) + i(c2(v^x+) + C!(v^x_)) + 

+^(F2(2a;2)-Fi(2xi)); (10) 

Q+ = n^df - 0^2) + Cihdi + C2hd2 + B; (ii) 

B = ^(c+(V2x+)C_(y2x„)+Fi(2xi) + F2(2x2)), (12) 

^It has been proved |12j that only for Laplacian (eUiptic) metric gik = Sik can the symmetry operators 
be reduced to second-order operators, and the corresponding Hamiltonians H^^\H^'^'' allow the separation 
of variables. 
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where the prime denotes the derivative of function with respect to its argument. A set of 
particular solutions of ^ was obtained in [T2], [HI, jU]- 

3. Two-dimensional models with twisted reducibility of 
supercharges 

In the previous Section it was shown that two different constructions with very different 
properties exist in two-dimensional SUSY QM. The first one (Subsection 2.1.) includes 
two scalar Hamiltonians H^^\ H^'^'^ (only one of them has normalizable^ zero-energy ground 
state wave function \E'o(a;) = exp(±^^)) and their 2x2 matrix partner hII\ The second 
one (Subsection 2.2.) contains only the scalar Hamiltonians H^'^\ H^^^ with no information 
about their ground-state energy in advance, and both H^^^ and H^^^ a priori (by construction) 
obey the important property of integrability with the symmetry operators R^^\ R^^^ of fourth 
order in derivatives (see ©)• The natural idea is to unite all the above tempting properties 
by combining these two constructions, i.e. by identifying the original Hamiltonian H^^^ as 
the same in both approaches. More precisely: let H^^^ of the form ^ have the superpartners 
H^l^ and if*^^-' in the first-order scheme, and at the same time the superpartner H^^'> in the 
second-order scheme. It is known ^2] that the simplest, reducible or quasi-factorizable, form 
of the second order supercharges = {Q^Y = qt'iii which is suitable for the construction 
described, leads to the i?— separation of variables, and therefore it is not considered here. 
All other models (excluding the case of elliptic metric gik = Sik in Q^) have been proved J2] 
not to be amenable to separation of variables; they have nontrivial fourth-order symmetry 
operators (jH)). The main idea to achieve the identification of H^^^ in the two approaches is 
to consider a class of models with second-order supercharges, which are quasi-factorizable, 
but with an intermediate twizst transformation (see also [21]): 

Q- = (g+)t = g+f/,,g-, (13) 
where Uik is a constant unitary matrix, q"^ were defined in and 

^For the case of supersymmetry not broken spontaneously. 
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with some new superpotential x- Such a generahzation of the notion of reducibihty (we 
shall call it twisted reducibility) is somehow reminiscent of the "gluing with shift" recipe in 
one-dimensional scalar 14| and matrix |2| HSUSY QM. The intertwining relations (jZj) with 
supercharges and the general expression for matrix Uik '■ 

U = aoC"o + i'^'a; + "c?^ = 1; a^, ai G R, 

(cTj are the Pauli matrices and (Jq is the unit matrix) give the system of four linear and one 



nonlinear equations for two functions X± = |(x i x) • 

ttsDX- + 2aidid2X- = 0; «inx+ - '2a3did2X+ = 0; (14) 

a2nx+ - 2q;o5i52X- = 0; aoDx- + 2a2did2X+ = 0; (15) 

{dkX-){dkX+) = 0, (16) 



where □ = 9f — Precisely the last equation (jlfij) is obviously most difficult to solve. Both 
the solutions of linear partial differential equations ()14|) -fll5 |) and the form of (fT^ depend 
crucially on the values of a, chosen. For the most sets of a's, including the general case with 
all aj 7^ as well as almost all degenerate cases with some vanishing, the corresponding 
potentials allow the separation of variables and are ignored here. 

The only exception to the above rule, and therefore the most interesting quantum models, 
corresponds to the case*^ when oq = ai = a2 = 0; a-^ 7^ 0, i.e. U = a^. Then, the metric 
of supercharges is hyperbolic, i.e. belong to the class discussed in Subsection 2.2. 
For these models (due to ()14|1 - (fTH|) only), the supercharges are represented in terms of four 
arbitrary real functions fii^2, '■ 

X- = /i+(x+) +/i_(x_), X+ = + f^2{x2), (17) 

Hence, the last equation (fT^ rewritten via = /i' takes the form of the functional equation: 

0i(xi)[0+(a:+) + 0_ (a;_)] + ^^2) [M^+) ~ <P-{x^)] = 0. (18) 

'^The system with ao = 0:2 = 0:3 = 0; ai 7^ 0, i.e. [/ = cti leads to analogous results with the 
substitution x± ^ Xt- 
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It is reasonable to formulate here the important specific property of solution ()17|). The 
superpotential 

Xix) =X+ + X-= l^iixi) + lJ,2ix2) + fi+{x+) + (19) 

leads to an expression for the quantum potential V^^\x) (see the first of Eq.®), which also 
has the form of the sum: 

l/(0)(f) = (^dix{x)y - hdfxix) = viixi) + V2{X2) + v+{x+) + t;_(x_); , (20) 

with t>i_2 = fJ'i\ — ^fJ'i2y "^i = fi'±—hfi'^. It may be seen that both terms in quantum potential 
()20|) separately have the form of the sums as in ()19p. Therefore, at the quasi-classical limit 
(see Sections 4 and 5 below), where only the first term (^dix{x)j survives, the potential 
is also represented in a form similar to (|20|) but with truncated vi^2,+ -- Both in the quantum 
and classical contexts, form ()2()|1 seems to be typical for a wide class of integrable two- 
dimensional models, considered within very different approaches in the literature (see j22j . 
|18j . as examples). This is why the following statement might be useful (at least, in the 
classical framework). Thus, if the general solution for the superpotential x{x) in relation 

Vd\^) = {dix{x)y = V,{xi) + V2{X2) + v+{x+) + v_{x^) (21) 

is of the form of (fTU|) . precisely the functional equation ()18|) must be fulfilled. This equation 
ensures the mutual cancellation of crossed terms in ()21|) and is therefore very important for 
this class of model. The general solution of ()18j] was found by D.Nishnianidze' (see 



- «0i,2 + ^01,2 + c; x = ± =; a, 6, c = const, (22) 

J ^a<p\ + b(p^ + c 

and explicit expressions for (l)± can be obtained from (fTH|l . Among the functions that satisfy 
conditions there exists a set of solutions of Eq.(fTH| possessing the periodicity property. 
For example, 

sn(ax|fc) cn{ax\k) dn{ax\k) 

(P+{x) = A — — ; (j)_{x) = A—- — — -— ; 23 

an[ax\K) /e^sn(aa;|fi;j cn(ax|/ej 

4>i{x) = 4'2{x) = B (hi{\/2ax\k) , 

where A,B,a are real constants and sn(-|A;), cn(-|fc) and dn(-|fc) are Jacobi elliptic functions 
with modulus k. They are doubly periodic on the complex plane of argument x, but 
'Private communication. 
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in the case k = 1 the real period becomes infinite and the elhptic functions turn into the 
hyperbohc functions sinh and cosh. Restricting ourselves in ()22j) to non-periodic functions 
on a whole plane {xi,X2) satisfying (j221), which do not happens in systems with separation 
of variables, two families of models exist (see [21] )• One of them is represented by the 
two-dimensional Morse potential, with (pi = (j)2 = Be~°'^; = 2 A; (/)_ = 2 A coth {ax/^/2): 



y(0) ^ (52g-2aa;i ^ fiBaC-'^''') + {B'^e 



■2ax2 



+2A{2A + h^) 



sinh 



hBae~ 
' + 8^2. 



0LX2 ^ 



(24) 



and the second by the two-dimensional Poschl- Teller potential with 



A [sinh (72 



ax 



52 



= B tanh {ax\. 

B{B+ha) 
cosh2(-^(xi+X2)) 

A-h\/2acosh {V2axi) 



B' 



B{B+ha) 
cosh^ {^{xi-X2)) 



+ 



+ 



(25) 



sinh^ (\/2axi) sinh^ {\/2ax2) 

Other members of these families can be obtained by using two discrete symmetries of solu- 
tions of Eq. iPIl : 

5*1: < 01 (Xi), 02(2:2), 0+(a;+),0-(x_) > y <^ 0+(Xl),0_(x2),0l(x+),02(X- 



S'2 : |0i(a;i), 02(0:2), 0+(x+),0_(2;_ 
and different combinations thereof. 



5 r + 



(26) 



4. Supersymmetric extensions of scalar Hamiltonians 
and their integrability 



In the previous Section we presented the explicit forms (j211), ()25|) of the terms in Eq.(j2ni) 
for a certain class of quantum integrable Hamiltonians. In this Section we shall build their 
classical and quantum SUSY extensions and we shall also demonstrate their integrability 
properties. 



4.1. The classical limit for H'^^^ 

First, we shall consider for H^^^ its classical limit H^^\ for which the integral of motion 
Rci exists: {H^^\r^^i}p = ({•, -jp denotes standard Poisson brackets). This can be done 



by the simple limit procedure ^ — > in Eq.Q. The practical recipe is as follows. One has 
to replace all operators —ihdi by momenta Pi and skip all derivatives of functions, which 
include h as a. multiplier. One thus obtains^: 

The integral of motion has the form R^^^ = Q^Q~^. Its involution with the Hamiltonian 
can be checked either by direct calculation or by a simpler two-step procedure, proposed in 
general form in ^j. It is instructive to perform it in the context of the models considered 
here. First, one has to prove the intermediate relations: 

{Hf,Q%}p = ±2z{<P'^ + <P'^)Q^,. (27) 

From the definition of Poisson brackets 

T 2^V2(0+(0i0;- 0202) + 0202)). (28) 

The last term can be transformed with the use of the functional equation ()18|1 : 

i=0_,9_ (0? + 08 + ^ 



+ \VmJl - (/>202) + 0-(010'l + 0202) = ^0+^- (01 + 08 + ^0-^+ (01 + 4>V) 



[9_(0+(0i + 02)2) + 9+(0_(0i - 02)2) - 20+9_(0i02) + 20_9+(0i02)] 
[-9_(0_(02 - 02)) - a+(0+(02 - 02)) - 20+9_(0i02) + 20_9+(0i02)] 

-^(0V + 0-)(0?-0^)- 



Substituting this into (j28|) . we prove (j27D. Then, ^ leads to the involution of R^^^ with 



(0) 

ItUlllig UlilO lliUU ljZ.Ojl ■ WC piUVC I 1/ . ±liCli, IjZ. I )l ICaUO UllC lliVUlUUlUli Ul 11, 

(0). 

(0) ^+^-1 _ ruW ^+1 ^- , r^+ru(0) 



It is easy to check that the same classical limit procedure for the second scalar Hamiltonian 
iJ'-^^ and for the components of the matrix H^l^ in Q leads to the simple results: H^f^ = 
JHere and below we will omit the arguments of (j)i^2 and 0±, which implies that 0i_2 = 01,2(2^1,2) and 
(f>± = (f>±{x±). 
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H^^^ and H^^^i = 6ikH^^\ Naturally, the corresponding integrals of motion coincide too: 
-^d'* ~ -^d"* ci ~ ^ikR^d ■ -'■^ Subsection we shall construct another classical 

limit of H that will also include Grassmanian dynamical variables in addition to pj and Xj, 
and this can be interpreted as a SUSY- extension of H^^\ 



4.2. The SUSY extension of classical scalar Hamiltonians 

It is well-known [3], jS] that the 2D representation of SUSY algebra, reviewed in Sec- 
tion 2, can be obtained by the canonical quantization from the classical system with the 
Hamiltonian: 

Hci = PiPi + {dix{x))idixix)) + 2tdAx{x)ijl^f, (29) 

where ip} and ^p'^ are Grassmanian anticommuting variables: {V'i")^'^} = {hj = l,---,d; 
a, (3 = 1,2). One can define the Poisson bracket on the phase space of the system with 
classical bosonic and fermionic variables |25l as follows: 



^ dpj dx^ dx^ dpj dtp" dipf 

such that the canonical brackets are {pj, Xj}p = 6ij; , ipj} p = iSij. Thus, the Hamiltonian 
is involved in the SUSY algebra: 

{Qa, Ql3}p = i^apHcl] {Hch Qa]p = 

with classical supercharges 

To quantize this model, one has to introduce the bosonic operators pi and Xj with canonical 
commutation relations, and fermionic ones t/'" obeying {tp"", ipj} = h5ij5"'^ . At this stage, 
it is convenient to introduce the fermionic operators 'ip'^ = (v^)~^(^j =F ii^j) with anticom- 
mutation relations {ip^,ipj} = h6ij and {4'^,ipj} = {V'l^^V'/} = 0. These can be treated 
as creation/annihilation operators in the system of d spin 1/2 fermions. Thus, the quantum 
counterpart of the Hamiltonian ()29|1 - the Superhamiltonian - in terms of these operators 
takes the form: 

H = -h%d, + [{djxMx) - md,x)] + 2{d,d,x)^Pti^J- (30) 

11 



Together with the quantum supercharges = ±i{V2)^'^{Q^ =1= iQ'^) = {±ipj + {djx))4'f 
this generates the algebra (H}. To reproduce the matrix form (0) - © one should choose 
the matrix representation for the creation and annihilation operators For d = 2, these 
operators are 4x4 matrices, and a possible choice is as follows [5 : ^'t = — -^4,3); 

ijj2 = /l^/^(-E3,i + -£^4,2) and = (V'/)"^ (here matrices Em,n are defined as {Em,n)ik = SmiSnj)- 
Thus, we have obtained a matrix realization of two-dimensional SUSY QM (Q) - © by means 
of the canonical quantization of a certain classical model ()29|) . 

One can see that in this representation H has a block-diagonal structure. The origin of 
this feature of the model is the conservation of the fermion number [H, N] = 0, with fermion 
number operator N = ^^ip~ . Therefore, each component of the Superhamiltonian acts in a 
space of states with a fixed fermion number. 

In our case d = 2, this structure is rather simple. Let us define a basis in the state space: 
the vacuum |00>, which is annihilated by , and the excited states |10>= h~^^'^ijjf\00>; 
|01>= h'^^'^tp2\00>; |11>= h~^^2''Pf\00> (the /i-dependent multipliers provide the proper 
normalization of the state vectors: <mn\mn>= 1, \fTn,n = 0,1). Thus H^^^ acts in one- 
dimensional space with a fermion number of and a single basis element |00>; for H^'^\ 
the fermion number is 2 and the basis element is |11>. The Hamiltonian acts in 
two-dimensional state space {|01>, |10>}, where the fermion number is 1. 

The conclusion to be drawn from this derivation is as follows. Having the classical system 
with H^^^ = pjPj + {djX){djx), one can construct its SUSY extension of the form (j^ . This 
classical SUSY extension can be quantized canonically to obtain the quantum Superhamil- 
tonian the original H^^^ being the classical limit of H^^^ - the first scalar component 
of the quantum Superhamiltonian. In our case, Hj!^^ was integrable (see previous Subsec- 
tion), and we shall explicitly find the integral of motion for its quantum SUSY-extension 
(Superhamiltonian). An analogous problem was investigated by alternative methods in |!26j . 
but for a much more narrow class of models (amenable to separation of variables). 

4.3. Integrals of motion for the quantum and classical SUSY ex- 
tensions 

We start from construction of the quantum integral of motion R for the Superhamiltonian 
()30|): [H,R] = 0, i.e. of conserved operators for each component of the Superhamilto- 
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nian: 



(31) 



= 0; (32) 

(note that the last commutator is the matrix one). The exphcit expression for R^^^ can be 
obtained from (jHJ and (fT^ : 

R'-'^ =qtU^kqk&mlqi. (33) 

One can obtain R^"^^ from R^^^ by the substitutions —qj and qf —qj (since 

turns to if*^^^ after the substitution x(x) ^ —x{^))'- 

R^'^ = qiU,,qtq-U^iq^. (34) 

With respect to i?*^^-* , the form of intertwining relations (0) tells us how to build this symmetry 
operator explicitly. One can check that the following matrix operator of sixth order in 
derivatives 

R^S = qiR^'\t. (35) 

satisfies Eq. ljH^ . 

It is clear from the material of the previous Subsection that knowledge of (i = 0, 1, 2) 
provides a symmetry operator R for the whole block- diagonal Superhamiltonian H. In order 
to construct an all-sector expression for R : 

R = RWpW + i?(2)p(2) + R^P^\ (36) 

we introduce the corresponding projection operators P^*^ (i = 0, 1,2) onto subspaces with 
definite fermion numbers. The scalar projectors P^'^^ and P^"^^ give unity when acting on 
|00> and |11>, respectively, and give zero otherwise. The components P^^^ {i,k = 1,2) of 
the 2x2 matrix projector operators P^^-* transform the k—th component of the state vector 
into its i—th component'^ and are zero on other components of the state vector. One can 
check directly that an explicit form of these operators P leads to: 

R = -R('^^h-^^tij^ijti^t-R^'^^~'i^ti^ti^ii^2-R^^^ 
-pgr^7/>rv>+7/>+v>2" + Pjyr + R^^^h-^ti^^, (37) 

'^By definition, "the first component" is |10>, and "the second component" is |01>. 
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where i?*^*-* are given by (jH^ - (jHS))- Naturally, ^7\i can be simplified by employing anti- 
commutation relations for tjjf. One should not be confused by the presence of the negative 
powers of the Planck constant in Eg. ()H7|1 since they disappear in all matrix elements for R. 

Let us prove straightforwardly that the R constructed commutes with H. The Super- 
hamiltonian can be presented similarly to (jHUj) : 

H = i7(0)p{0) + ^(2)p(2) ^ ifWpW^ 

where H^'^ are given by Eqs.©. By definition [p(o),p(2)] = [P/;J\ P^^)] ^ [P^^^ p(o)] = 0, 
and therefore, due to (p?T|l . we have: 

iff./?i = i4'/^f,i?«p<;']. 

Employing the explicit form ()37j) of P*-*^ and Eq. ()H2|l ). one can see that the commutator in 
the rhs vanishes, completing our proof. 

From expression (jSTj), its classical limit Rd can be constructed by means of the substi- 
tution h~^^'^ijj°' — > -ipj" for fermionic operators. Finally, the total classical integral of motion 
reads: 

Rci = + '^{R^'L - R'iLdmi - i'li'l) + \{R'iL - RtLdmi + ^l^D + 
+^(^!2cz + R'ilMi'l + ^iV2) + {R'lli + R'iL - 2Rf)i^\i^i^ii^i (38) 

where its components are: 

Rd = i?I? = (p?-p^ + 0?-0^-20+0_)' + 2((0+ + 0_)pi-(0+-0_)p2)'; 

(0). 

d ' 



R^Ki = (p? + 0? + ^(0U0'_)+2v^0+0_^P 
R^li = (ipi + h + ^(0+ + 0-)) (^-W2 + 02 + ^(0+ - 0-) ) R^^^; 



(39) 



4ld = ( *P2 + 02 + ;^(0+ - 0-)) (-^Pi + 01 + ^(0+ + 0-) ) 
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5. The flipped potentials and the classical Hamilton- 
Jacobi equation 



In this Section we shall establish links between the Hamilton- Jacobi equations of Classical 
Mechanics and the equation for the superpotential. Starting from Eq. ()21|1 . one can see that 
the condition for the classical Hamiltonian 

Hci = Pi + V{x) 

to be supersymmetric with superpotential x{^) takes the form: 

dxix) dxix) 



V{x) 



(40) 



dxi dxi 

On the other hand, for the Hamiltonian h^i = pf + U (x) with potential U (x) and the classical 
action functional S, the well known Hamilton- Jacobi equation j22j reads: 

dS 



OS ^ dS dS 
dt ^' dxi ' dx2 



dxd 



Xi,X2, 



, Xd) = 0. 



(41) 



There being no explicit dependence on time in hd, one looks for its solutions of the form 
S(t,xi,X2, - ■ ■ ,Xd) = W{xi,X2, ■ ■ ■ ,Xd) — Et, and the time-independent Hamilton- Jacobi 
equation becomes: 

E = (diW)'^ + U{x), (42) 

where W{xi, X2, - ■ ■ , Xd) is usually referred to as the Hamilton characteristic function. Solu- 
tions of (jl^ in the case E = are obviously connected with those of (pn|) : 



X{x) = ±iW{x). 



(43) 



Eq. ()42|) . with zero energy, can alternatively be thought of as a condition for the "flipped" 
classical potential V = —U to be supersymmetric, i.e. V should satisfy (j^ . with x and W 
related by Thus, we find 

W{xi,X2) = + Ai2(x2) + /i+(a;+) + /U-(a;_)] 

as the Hamilton characteristic function of the system. The system of equations of motion 
for E = 



Xi 



X2 



dW 
dxi 
dW 

dX2 



1 
1 

7!^ 



+ 



1 
1 

T2' 



(44) 
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is not amenable to separation of variables and in general has (non-physical) complex solu- 
tions. One may see that system ()44|) becomes real, and bona fide solutions exist for the 
specific complexification of the Poschl- Teller model ()25j) : namely, with purely imaginary a. 
In contrast to the complexification of two-dimensional Morse potential in ||2£,, this one is 
PT- invariant. 

Analogous classical systems with "fiipped" potentials were investigated in for the 
case of (i = 2 integrable models of the Liouville type. For these systems the Hamilton- Jacobi 
equations were separable in elliptic, polar, parabolic and Cartesian coordinates. The struc- 
ture of related supersymmetric models (also with separation of variables) in the quantum 
domain has been investigated in ^H] via canonical quantization. In particular, it was shown 
that there are two essentially different supersymmetric extensions (two different superpoten- 
tials) for a given separable classical solution of the Hamilton- Jacobi equation. In the present 
paper our strategy is just the opposite. Namely, starting from scalar quantum d = 2 sys- 
tems which do not allow for separation of variables, but do have non-trivial symmetry 
operators, we construct their quantum SUSY-extension. Then, we describe corresponding 
classical SUSY-ext ended systems and their integrals of motion. Finally, the link between 
this kind of classical system and the Hamilton- Jacobi approach for related systems with 
"flipped" potential provides the integrability of these "flipped" systems too. The necessary 
explicit expressions for integrals of motion can be obtained from by the substi- 

tution (f) —>■ ±i<f), which is equivalent to PHj) . It should be remarked that (analogously to 
[T^ ) besides an arbitrary common sign in (pSj) there is the additional non-uniqueness of the 
superpotential for this class of model. Indeed, equation pUjl has two independent solutions 
for flxed original classical potential: x{^) ^ind x{^)- To prove this statement, one has to 
check that (9;x)^ = {dixY X = /^i + + Ai+ + A*- and X = l^i + 1^2 — 1^+ — 1^- (see (fTTj) ) 
due to the nonlinear equation |T6|l . 
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